Fokker-Planck equation approach to vehicle statistics 
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This contribution presents a derivation of the steady-state distribution of velocities and distances 
of vehicles in freeway traffic which has been suggested for the evaluation of interaction potentials 
among vehicles (see preprint cond- mat/0301484 1. Despite the forwardly directed interactions and the 
additional driving terms in vehicle traffic, the steady-state velocity and distance distributions agree 
with the equilibrium distributions of classical many-particle systems with symmetrical interactions, 
if the system is large enough. Finally, this analytical result is confirmed by computer simulations. 

In the particular driven-many particle system we discuss, driver-vehicle units play the role of the particles. Here, 
we will describe their behavior by the coupled car-following equations 



dvj _ vp - Vj 
dt ~ t 



+ /(*<) -7/(*i-i)+&(*) 



(1) 



where Vi(t) = dfi/dt is the speed of vehicle i at time t, vq the maximum velocity, s,(t) = ri(t) — n+iit) the distance, 
and represents a white noise fluctuation term. The term 7/(s,_x) with < 7 < 1 allows to study different cases: 
7 = corresponds to the case of forwardly directed interactions of vehicles, while 7=1 corresponds to symmetrical 
interactions of classical particles in forward and backward direction fulfilling the physical law of "actio = reactio" . 

The above stochastic differential equation (Langevin equation) can be rewritten in terms of an equivalent Fokker- 
Planck equation. With the definitions 

W( Si ) = « +r[/( Si )-7/(si-i)], 
dU( Si ) 



f(Si) 



ds. 



&(*)> = 0, 
&(*)&(*0> = D5 l3 5{t-t'), 

and P = P(si, s n ,vi, ...,v n ,t) , 



this Fokker-Planck equation reads 
OF 



9 \l \ Dl 9 

—dsi I dt 



( W( Si )-Vi 



P 



Dd 2 P 



(2) 



(3) 



where we assume periodic boundary conditions Vk+ n (t) — Vk(t) and Sk+ n {t) = Sk(t) for a freeway of length L. In the 
following, we will show that the ansatz 



P( Sl , ...,s n ,v 1 ,...,v n )= Afc-^A^ye+Bs^-E^-vf/m 



(4) 



is a stationary solution of the above Fokker-Planck equation, if the parameters V and 9 are properly chosen. The 
parameter B is required to specify the actual vehicle density (i.e. to ensure J^j s » = 
In Eq. @, 



N= ids,... I ds n I dm... I dv n e-ZA u ^y e+B ^e-Z^-vr/W) 
is the normalization constant, 

V(t) = (Vi) = I dsx . . . I ds n I dvx . . . I dv n ViP(si, . . . , s n ,vi, . . .,v n ,t) 



(5) 



(6) 
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is the average vehicle velocity, and 



6{t) = ({n l -V?) 



dsi 



ds n I dv% 



dv n ( Vi - V) 2 P( Sl 



(7) 



the velocity variance. In the following, we will restrict our investigation to the stationary case with dV/dt = and 
d6/dt = 0, which presupposes that the instability condition of Eq. Q is not fulfilled. For traffic systems with 7 = 0, 
the instability condition is known to be of the form 



dW(s) 1 
ds > 2r 



(8) 



If this condition applies, stop-and-go traffic will emerge. 
Differentiation of Q gives: 



d 



E 7)7}( Vi ~ = E(^ - 



1 dU{ Sl ) 
6 ds,, 



B 



P = E^' 1 - v i+l) 



B- 



f(Si 



p. 



(9) 



E 



d_ fW( Si )-Vi 
dv,. 



P 



p = Et+E 



v t - W( Si ) 



(«< - V) 



P 



and 



D d 2 P _D ^ 



vt - V 



2 dv,, 2 2 

We will now insert this into Eq. and use the fact that 

E^±i p = E^ 



p. 



(10) 



(11) 



(12) 



for any i-dependent variable gi, i.e. indices can be shifted because of the assumed periodic boundary conditions. In 
this way we find 



dP 
~dt 



DP 
~26 



-E 



Vo - Vi 



+ /(*)- 7/(*i-i) 



26» 2 



Ansatz Q can only be a stationary solution with dP/dt = 0, if 

1 _ 2 
6 ~ ~Dt ' 



(13) 



(14) 



which relates to the fluctuation-dissipation theorem. With this, (v;+i — vi) = (v;+i — V) — («, — V), and (vq — Vi) 
(vo-V)-(Vi-V)],we&ad 



dP 



at 



We will distinguish the following cases: 



l^(vo-V)( Vl -V) 



P. 



(15) 



1. In the case of a classical many-particle system with momentum conservation (7 = 1) and energy conservation, 
i.e. no driving (vo — 0) and no dissipation (r — + 00), we find dP/dt = 0, i.e. ansatz Q is an exact stationary 
solution of the Fokker-Planck equation @ . 



2. In the case of vehicle traffic (7 = 0), we have to show that the additional term 



ij>n-i-v) 



f(Si) + 



v -V 



(16) 



3 



disappears (where we have again shifted indices). Let us first note that, with the factorization assumption, we 
can state 



lim — V7ui+i - V) 

n — >rv~i 71 * * 



/(*) 



vo-V 



lim < lim — 

n— >oc Tl A — ' n— >oo fi — * 



f{Si 



v -V 



(17) 



The first factor vanishes because of V = linin^oo ~Yli v ii but ^ ne second factor disappears as well: Dividing 
Eq. P by n and summing up over i gives 



1 x - dvi 1 x 
n^—'dt ~ n 2—> 



v - Vi 1 



(18) 



In the limit n — ► oo of large enough particle numbers n, the left-hand side converges to dV/dt, while the 
last term on the right-hand side converges to 0. In the assumed stationary case with dV/dt = and using 
vo — Vi — (vq — V) — fa — V), this implies 



lim — 



vo - v t 



+ /fa) 



= and lim — 

tj. — >r>Ci 11 ' ^ 



v - V 



+ f(Si) 



= 0. 



(19) 



Nevertheless, it is not obvious what happens for finite systems, as the standard deviation of 53ifa+i — ^0 ^ s Vnd- In 
order to see how the single-particle distributions depend on n, let us again assume a factorizing solution 



P(sx, ■ ■ -,a n ,v 1 , . . -,v n ,t) = Y\g{s l ,t) Y\ h(vj,t) 



i=l 



in generalization of ansatz Inserting this into l|15|l with 7 = gives 

9P _sr TtlaMHvi+ut)] _ 1 Vr _ . 



or 



/(*) + 



/fa) + 



P 



g(.3i,t)h(v l+ i,t) . 



(20) 



(21) 



(22) 



The left-hand side represents the average temporal change of the one-particle distribution functions 3 fa, t) of the netto 
distance and /ifa+i, t) of the speed, while the right-hand side converges to zero with growing system size according to 
the central limit theorem and the factorization assumed with Eq. 120[) . i.e. the statistical independence of the variables 
Vi and Si. That is, while is an exact equilibrium solution for a classical many-particle system, it is also expected to 
be a steady-state solution of driven many-particle systems of the kind QJ , even if the potential is forwardly directed 
rather than symmetric! 

Note that, because of the factorization ansatz (@J and pU)l. one may use the approximation 



l£(„ <+1 _V) 

71 ^ ^ 



/fa 



vo-V 



-Y,(v i+1 -v) j i E 



/fa 



vo-V 



(23) 



where both factors on the right-hand side converge to zero because of l|19fl . Therefore, the convergence should be 
particularly fast. Moreover, in empirical evaluations, the estimator of V is i- e - t ne ^ rs ^ factor on the right- 

hand becomes exactly zero. For all these reasons, it is expected that 



g(s) oc e -^/ e+Bs l 
is a good approximation of the empirical distance distribution and 

h(v) cx e-^ 2 ^ 



(24) 



(25) 



a good approximation of the empirical velocity distribution, where V = — Vi, 8 = — Y)f fa — V) 2 , and n > 50. 
This is actually confirmed by numerical simulations of Eq. see Figs. ^ and [21 The numerical results (symbols) 
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FIG. 1: Semilogarithmic plot of the frequency distributions of the velocity (left) and the distance (right) according to numerical 
simulations for symmetrical interactions (top) and forwardly-directed ones (below). The assumed interaction potential was 
U(s) oc s _1 for s > and otherwise. 



agree with the steady-state solution (solid line). Within the statistical variation and apart from finite size corrections, 
the results are the same for symmetrical and forwardly directed interaction potentials. 
Finally, let us investigate the Hamiltonian 

W = T + V = £fc^ + ]Tt/( Sl ). (26) 

i i 

If dV/dt ~ 0, we can derive the following relations: 

dH _ dT dV 
dt dt dt 

dvi \r^dll(si) ( dsi dri dsi drj+i 



dt dsi \dri dt dri + \ dt 

ivi 
~d7 



- v) + f( Si ) - 7 /(*i-i) + m) + E - «i) 

i ^ i 

i i 

Y^i^izlL + ^-VMt)- (27) 
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FIG. 2: Semilogarithmic plot of the frequency distributions of the velocity (left) and the distance (right) according to numerical 
simulations for symmetrical interactions (top) and forwardly-directed ones (below). The applied interaction potential was 
U(s) oc s~ 4 for s > and otherwise. 



Comparing this with (|15fl shows that 

dP 

~dt ' 9 dt 



dP PdH ly, 



(iH - V f 



-(Vi-V)^(t) 



P. 



Correspondingly, we have 

dU 
~dt 



E 



(Vi - Vf 



Vi-V 



(28) 



(29) 



in the stationary state dP/dt = 0. We will again distinguish two different cases: 

1. In a conservative system with no fluctuations (£i(t) = = D) and no dissipation (r — ► oo), we have dTL/dt = 0, 
independently of whether the interactions are symmetric or forwardly directed. 

2. For many-particle systems with fluctuation terms and/or dissipation, one can show 

' 1 d(vi - Vf ' 



Uvt-V)) = (± d{Vl dt V) )- VJL ^Y.^ - V ) + 7<fo - V ) 2 ) - - 7/K-i)1K - V)) 



2 dt 

IdO v -V 
2~di t 



-V) + -- (f( Si ) - 7/(«i-i)>««i> - V) . 



(30) 



This can be found by multiplication of Eq. (JIJ with (v{ — V) and calculation of the ensemble average, using 
the factorization ansatz JIJ or (|20|) . The first term on the right-hand side vanishes under the assumption of a 



6 



stationary state. The second and the fourth term vanish because of (Uj) = V. Therefore, 



and 



1 



vf 



(31) 



r 



if n is large enough. Without dissipation (r — > oo), (£i(«i — V)) becomes zero, while it is finite otherwise. 
Together with Eq. I)2H[I. we arrive at 



i.e. in the statistical average we have dH/dt — 0. This is also expected for systems with many particles. 

In conclusion, the equilibrium solution of conservative many-particle systems is a good approximation for steady- 
state solutions (dP/dt) of driven many-particle systems of the kind with asymmetrical interactions, driving and 
dissipation effects, if the system is large enough, i.e. n>l. For small systems, we expect that fluctuations become 
essential. A more detailed elaboration is presently in preparation and will be submitted, soon. 




(32) 



